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ABSTRACT 

Context. In 2007 we reported two cometary shaped sources in the vicinity of Sgr A* (0.8" and 3.4" projected distance), 
named X7 and X3. The symmetry axes of the two sources are aligned to within 5° in the plane of the sky and the tips 
of their bow-shocks point towards Sgr A*. Our measurements show that the proper motion vectors of both features are 
pointing in directions more than 45° away from the line that connects them with Sgr A*. This misalignment of the 
bow-shock symmetry axes and their proper motion vectors, combined with the high proper motion velocities of several 
lOOkms"^, suggest that the bow-shocks must be produced by an interaction with some external fast wind, possibly 
coming from Sgr A*, or stars in its vicinity. 

Aims. We have developed a bow-shock model to fit the observed morphology and constrain the source of the external 
wind. 

Methods. The result of our modeling gives the best solution for bow-shock standoff distances for the two features, which 
allows us to estimate the velocity of the external wind, making sure that all likely stellar types of the bow-shock stars 
are considered. 

Results. We show that neither of the two bow-shocks (one of which is clearly associated with a stellar source) can be 
produced by influence of a stellar wind of a single mass-losing star in the central parsec. Instead, an outflow carrying a 
momentum comparable to the one contributed by the ensemble of the massive young stars, can drive shock velocities 
capable of producing the observed cometary features. We argue that a coUimated outflow arising perpendicular to the 
plane of the clockwise rotating stars (CWS), can easily account for the two features and the mini-cavity. However, the 
collective wind from the CWS has a scale of > 10". The presence of a strong, mass-loaded outbound wind at projected 
distances from Sgr A* of < 1" is in fact in agreement with models that predict a highly inefficient accretion onto the 
central black hole due to a strongly radius dependent accretion flow. 
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1. Introduction 

Analyses of stellar orbits in the central arcsecond of 
the Milky Way have provided indisputable evidence that 
the central object at the position of the radio source 



(SMBH; e.g. 


Eckart et al. 


20021 ISchodel et al. 


2003 


Eisenhauer et ai.||2005| |Ghez et al.||2005l |2008t Gil 


iessenl 


et al. 2009 


). Sgr A* is located at the distance of ~8kpc, 



oi - 4x WMq. 
There is long-standing evidence for the existence of out- 
flow(s) at the Galactic Center (GC). The central parsec 
of the Milky W ay harbors a cluster of massive stars (e.g. 



Paumard et al. 2006J that contribut es ~ 3 x 10 Mp^ yr 



to the center in form of stellar winds ( [Najarro et al.|1997 ). 
However, less than 1% of this material is actually accreted 



onto the SMBH (Baganoff et al. 2003; Bower et al. 2003 



Marrone et al.|2006( ) and thus has to be expelled out of the 
center. 



* Present address: Department of Astronomy and 
Astrophysics, University of Toronto, 50 St. George Street, 
Toronto ON M5S 3H4, Canada 



The streamers of gas and dust in the central few par- 
sees of the Galaxy show a bubble-like region of lower density 
(called mini-cavity), located ~3.5" southwest of Sgr A*. It 
was first po inted out on cm-radio maps by |Yusef-Zadeh 



et al. (1990). The strong Fe[ni] line emission seen toward 
that region fEckart et al. 1992 Lutz et al. 1993) is con- 



sistent with gas excited by a collision with a fast (>1000 
kms~^) wind from a source within the central few arcsec- 



1992) 



onds (Yusef-Zadeh & Wardle 1993 lYusef-Zadeh & Meha 



In |Muzic et aL] ( |2007[ ) we presented NACO/VLT multi- 
epoch oEservatlons at 3.8 fim (L'-band) which allowed us to 
derive proper motions of narrow filamentary features asso- 
ciated with the gas and dust streamers of the mini-spiral. 
Proper motions of several bow-shock sources have also been 
reported. The analysis has shown that the shape and the 
motion of the filaments is not consistent with a purely 
Keplerian motion of the gas in the potential of Sgr A* and 
that additional mechanisms must be responsible for their 
formation and motion. We argued that the properties of 
the filaments are probably related to an outflow from the 
disk of young mass- losing stars around Sgr A*. In part, the 
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Fig. 1. NACO L'-band (3.8 /im) image of the Galactic Center region showing projected positions of SgrA* and the 
cometary-shaped features X3 and X7. The arrows show proper motions. North is up, east is to the left. 



outflow may originate from the immediate vicinity of the 
black hole itself. 



In Muzic et al. 
cometary-shaped 



( 2007 1 we report the existence of the two 
features X3 and XT, located at projected 
distances of 3.4" and 0.8" from SgrA*, respectively. The 
symmetry axes of the two bow shocks are almost aligned 
(within 5°) and point towards SgrA* (see Fig.jl]). At the 
same time, their proper motion vectors are not pointing 
in the direction of the symmetry axes, as it would be ex- 
pected if the bow-shock shape were produced by a super- 
sonic motion of a mass-losing star through a static interstel- 
lar medium. In general, bow shock appearance can also be 
produced by an external supersonic wind interacting with 
a slower wind from a mass-losing star. If the velocity con- 
trast of the two winds is large, one will observe a bow shock 
pointing towards the source of the external wind. If in ad- 
dition the star moves in a direction different than the po- 
sition of the external wind, the bow shock will point in the 
direction of the relative velocity between the star and the 
incident medium. 

In this paper we present new proper motion measure- 
ments for the two cometary-shaped sources X3 and X7 
(section[3|. In section[4]we explain details of the bow-shock 
modeling that was used to fit those features (sectionlsl) . 
Possible thee-dimensional positions of the sources are dis- 
cussed in section[6] In section [7] we argue about the nature 
of the two sources and the external wind held responsible 
for the observed arrangement. Summary and conclusions 
are given in section[8] 



2. Observations and Data Reduction 

The L' (3.8 fim) images were taken with 
NAOS/CONICA adaptive optics assisted 
ager/ spectrometer (Lenzen et al 



1998 



1998[ prandner et al.,206^[ ) at the UT 



Rousset 
( Yepun) 



et al. 

at the 



ESO VLT. The data set includes images from 7 epochs 
(2002.66, 2003.36, 2004.32, 2005.36, 2006.41, 2007.25 and 
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Fig. 2. L'-band proper motions of the two cometary fea- 
tures X3 and X7. The error bars show the la uncertainty 
of each measurement. The x and y components refer to RA 
and Dec, respectively. 

2008.40) with a resolution of ~100mas and a pixel scale 
of 27 mas/pixel. Data reduction (bad pixel correction, sky 
subtraction, flat field correction) and formation of final 
mosaics was performed using the DPUSER so ftware for 



astronomical image analysis (T. Ott; see also Eckart & 
Duhoux|p90| . 



3. Proper motions 

In order to obtain proper motions of extended features we 
have used t he same method described in detail in^luli31 

The 



the et al. 



( 2007 ) . Our measurements are shown in Fig. [2] 

distance to the GC is assumed to be 8 kpc. The results ob- 
tained from the L'-band data between 2002 and 2008, are as 
follows. X7: v„ = (-52±12)kms-\ V5=(546 ± 15) kms'^. 
X3: Va = (-95±13) kms"\ V5=(122± 14) kms'^ The un- 
certainties are 1 a uncertainties of the weighted linear fit to 
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Fig. 3. The two-dimensional b ow-shoc k shap e. Full line 
shows the analytic model from IWi lkinI ([19961 see Eq.[T| . 



Dotted lines show narrow solutions ( Zhang &: Zheng|l99T I 
for two coUimation angles 9q—A0° and 610=60" . The axes 
are in units of Rq; the star is at the center of the coordi- 
nate system. 

Inset: 3D bow-shock sketch, explaining the conventions for 
the inclination angle used in this work. i=0° when the bow- 
shock symmetry axis is parallel to the line of sight (i.e 
we are looking directly at the apex). i=90° if the bow- 
shock symmetry axis is perpendicular to the line of sight. 
For a bow-shock inclined away from the observer (the ob- 
server looks at the envelope, but cannot see the apex), 
90° < i <180°. 



the positions vs. time. The proper motion velocity vectors 
of both sources are oriented in the northwest quadrant (see 
Fig. IT]) . We note here that t he results are i n agre ement with 



of IMuzic et al. 



where 



J 



20071, but there 



otted as moving 



results given in Tablet, 
is an error in their Fig. 5. 
towards north-east. 

The L-band feature X7 is coincident (in proje ction) 



2009). The 



with a K-band source S 50 ( Gillessen et al. 
reported proper motion of S 50 is Vq,=(-108 ± 5) kms~^, 
and V5=(361 ± 7) kms^^. While the orientations of the two 
proper motions are in agreement, there is a significant dif- 
ference in proper motion magnitudes. This, however, might 
give us an idea about the systematic uncertainties of the 
proper motions derived from the L-band extended sources. 
It seems reasonable to argue that we are not dealing with 
an accidental superposition of a stellar source with a dust 
blob along the line of sight, but that the L-band feature is 
indeed associated with the stellar source at the same posi- 
tion. 



4. Model 

4.1. Bow-shock shape 
4.1.1. Analytic model 



Wilkin ( 1996 ) derived a fully analytic solution for the shape 
of a bow shock produced by a star moving through the 
interstellar medium at supersonic velocity (see Fig. pi) . The 
two-dimensional shell shape is given by: 



R{e) = Ro CSC 9^/3(1 - 6* cot ( 



(1) 



where 9 is the polar angle from the axis of symmetry, as 
seen by the star at the coordinate origin. Rq is the so-called 
standoff distance obtained by balancing the ram pressures 



of the stellar wind and ambient medium at 0—0 and is given 
by: 



Rq 



(2) 



Here rhw is the stellar mass-loss rate, Vw the terminal veloc- 
ity of the stellar wind, Pa the ambient medium mass density, 
and Va is the velocity at which the star moves through the 
medium (i.e. the relative velocity of the ambient medium 
and the star velocity in case that the ambient medium is 
not stationary). 

We a ssume that the shell ha s a thickness as shown in 
Fig 3b in Mac Low et al. (1991). Then we rotate the two- 
dimensional shape around its axis of symmetry in order to 
obtain the three-dimensional shell. 



4.1.2. Narrow solution 



The model by Wilkin ( 1996 1 incorporates the instantaneous 
cooling approximation: the interaction between the stellar 
wind and the ambient medium takes place in an infinitely 
thin layer, the two fiows are fully mixed and immediately 
cooled. In this case Rq directly gives the distance of the 
star to the a pex of the bow shock. However, |Comer6n fc| 



Kaper 
of the 



( |1998j ) note that this might not be true if cooling 
shocked stellar wind is inefficient. In this case it is 
expected that the bow shock wo uld be located at a dis- 
tance somewhat greater than Rq (Comeron fc Kaper||1998 



Raga et al.|1997l|Povich et al.|2008p .TOT some GC sources , 
Tanner et al.| ( 2005 ) indeed reported that a better fit can be 
obtained if the apex of the bow shock is shifted away from 
the star. We have observed the same behavior for several 
other Galactic Center bow shocks. 



Zhang & Zheng ( 1997 ) have investigated the case where 
the wind ejection from the star is not necessarily isotropic; 
it is rather confined in a cone of solid angle — 27r(l — 
cos^o); where 29o is the opening angle in which the matter 
is ejected. The standoff distance is then given by: 



Ro 



(3) 



As can be seen in Fig.[3| with 9q getting smaller, the bow 
shock profile will have a narrower shape. For 6*0 — > 37r/4 the 
shape is very similar to the one given by Eq.[lj Although 
this model has been developed for the case of a collimated 
stellar wind, we do not assume such a coUimation, since it 
would imply that the coUimation of the wind is in the direc- 
tion of motion of the star, relative to the ambient medium. 
There is no a priori reason for the stellar velocity and wind 
coUimation direction to be linked. We introduce 6*0 as an 
additional factor that allows one to control the width of the 
bow shock outline, while causing the apex of the bow shock 
to be slightly displaced away from the star. 

4.2. Generating emission maps 

To generate a simulated observation, we start by illuminat- 
ing the shell by the star placed at the origin of the coor- 
dinate system and then calculate the emission along each 
ray intersecting the shell. The shell is inclined to the line of 
sight by the angle i, and also rotated by a position angle PA 
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in the plane parallel to the plane of the sky. Inset in Fig. [3] 
explains the convention used for the inclination angle. PA 
is measured East of North. 

Each parcel of the shell is assigned the optical depth r, 
calculated as: 

t(A) = TabsW + TscaW = ^ nd{a)Cext{a, X)da. (4) 

We assume a graphite+silicat e mixture with a po wer-law 
grain size distribution nj,{a) (Mathis et al. 1977), where 
a is the grain size. Dust extinction coefficient is Cext = 
T^a^iQabs + Qsca), where Qabs and Qsca ar e dust absorp 



tion and scattering efficiencies, respectively (Laor & Draine 
1993| L is the length of the shell parcel along the line of 
sight. 

The source X7 is polarized (see 
suggesting that scattering of the stef 



Muzic et aL||2007[ ), 
ar emission by dust 



particles in the bow-shock envelope is probably impor- 
tant. In order to account for scattering, we proceed 
in the following way: Emission of each parcel of the 
shell has contributions both from scattering and thermal 
emission: L=Lsca+Lth- Lm (xB{Td){l - e"''»''=)et/j and 
Lsca ofi~^€scaP{Qsca)s.~^""' , whcre d is the distance from 
the star, B{Td) is the black body emission at the dust tem- 
perature, integrated over grain sizes, and over wavelengths 
in our observing band, eth and Csca are thermal emission 
and scattering efficiencies, respectively. Dust temperature 
at a distance d (in parsecs) fro m the star can be calculated 

-1/6 



as Td = 27 On 



' lI% d~c'^ K (|Van Buren & McCray 



1988 



Kriigel, R]|2003 ), where a^m is dust grain size in ^m, anc 



L*^38 is stellar luminosity in lO^^ergss"^. P{9sca) is the 
(normalized) scattering function that controls the amount 
of forward scattering, and is given by: 



P{esca) 



1 + .g2 - 2gcos{0sca) ' 



(5) 



For (7 = 0, scattering is isotropic and P does not depend 
on the scattering angle 9 sea- 5 = 1 means full forward scat- 
tering. We have varied g between and 0.8, in steps of 
0.2. Since eth and eg^a are in general not known, we have 
chosen to test three different cases: {eth, Csca) = (1-0, 0.0), 
(0.0, 1.0), and (0.5, 0.5), i.e. thermal emission only, scat- 
tering only, and equal contributions from both processes. 
The shape of the thermal emission depends on the dust 
temperature distribution: we investigate central stars with 
luminosities (10^, 10^, 10'*)Lq. 

Finally, the resulting projection of the 3D-geometry 
onto the plane of the sky is rebinned to the pixel-scale of the 
NACO images and smoothed with a gaussian PSF having 
FWHM equivalent to the angular resolution of our images. 



5. Results 

5.1. X7 

Fig.|4] shows the best-fit results of our bow-shock modeling 
for the feature X7. We have tested the three different com- 
binations of {eth, ^sca), and five values of g={0.0, 0.2, 0.4, 
0.6, 0.8). In Fig.llj however, we show only those values that 
result in good fits, except for the left-most panels (blue 
contours), which we show to illustrate the behaviour of 



purely thermal emission. Thermal emission is dominant in 
the vicinity of the star, but cannot fit the extended tail that 
we observe in X7, unless the central star is much brighter 
than the cases investigated here. This, however, is not plau- 
sible since X7 is very faint in the K-band (see discussion 
below). The tail of the bow-shock is much better described 
by scattering. For {eth, esco) = (0.5, 0.5), we show only the 
solution for L* = 10^Lo. We choose to do so because the dif- 
ference in resulting contours for the three stellar luminosi- 
ties is small, and gives the same solution for the physically 
most interesting parameter, Rq. In cases when scattering 
is important (csco >0.5), more forward scattering (large 
g) results in a more compact model, thus not fitting well 
the outer contours. There is no significant influence of the 
parameters on the inner contours, due to relatively small 
size of the feature, and smoothing. These differences can be 
better observed in case of X3. 

For each set of parameters {eth, ^sca, g), we have tested 
different values of Rq. By changing Rq, the model preserves 
the same shape of the contours, but is as a whole expanded 
or shrunken. Therefore in Fig. [4] we plot two values of Rq 
for each set of parameters. We chose values in the way that 
shows how the change in Rq affects the fit. By changing its 
value by larger amount, the fit becomes inadequate. 

The best solutions for different sets of parameters are 
obtained for Rq w 2.5-lO^^cm, and PA = 50°, measured east 
of north. Note that in this case the best results are obtained 
using the simple analytic two-dimensional solution (eq.[T]). 
Both X7 and X3 have an unusually narrow appearance, and 
therefore are best fitted with inclination angles close to 90° . 

X7 coincides with a point source at shorter wave- 
lengths (see discussion on proper motions in sec- 
tion [3]). Phot ometric measuremen ts give H=18.9±0.1 and 
K=16^9±0.1 dSchodel et "al]|2010[ ). For the l ocal extinction 
at th e position of X7 we assume Ax=2.5 (Schodel et al. 
20T0|. 



In section [7T] we discuss possible stellar types and im- 
plications this has on the external wind parameters. 

5.2. X3 

Fig. [5] shows the best-fit results of our bow-shock model- 
ing for the feature X3. This feature is very elongated and a 
satisfactory fit cannot be obtained using the an alytic 2D so- 
lution. It requires a narrow model (see Section 4.1.2), with 
small opening angles Oq. As in the case of X7, the outer con- 
tours are represented better in models with lower g, while 
larger g values result in more compact inner contours. Here 
it is even more evident that thermal emission gives a too 
compact model. The elongated tail of X3 can only be well 
fitted with models that include scattering. Therefore we 
show only these solutions, in pairs of two different values 
of i?o- The best fit solutions give Rq «1.5-10^^cm, with 
6»o=30°, i = 90°, and PA = 55°. 

In contrast to X7, there is no detectable point source at 
the position of X3 in our Ks-band images. Local extinc tion 



at the position of X3 is Ak w2.7 (Schodel et al. 2010). In 



section 7.2 we discuss the possible nature of this source and 



^ available at http://www.astro.princeton.edu/~draine/ 



the implications this has on the external wind parameters. 



6. 3-dimensional arrangement 

Fig. [6] shows a 3D reconstruction of some of the fea- 
tures found in the central parsec of the Galaxy. The 
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Fig. 4. Results of our modehng for the feature X7 (black contours), observed in the epoch 2003.36 (the image is previously 
deconvolved and beam-restored to the nominal resolution of our L'-band images). Color contours represent the bow-shock 
model projected onto the plane of the sky. PA = 50° (E of N), and the shape is given by the analytic solution (section 
4.1.1 ). We show purely thermal emission in blue (left), pure scattering in orange (middle), and equal contributions from 
both processes in purple (right). For each set of parameters, we show the results for two different values of i?Oi in order to 
depict how change in the bow-shock size affects the fit. Contour levels are at (20, 30, 40, 50, 70, 90)% of the maximum. 
North is up, east is to the left. The axes are labeled in pixel offsets from the brightest point in the image; 1 pixel is 
equivalent to 0.027 mas. 
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Fig. 5. Best fit results of the modeling for the feature X3 (black contours), observed in the epoch 2003.36 (the image is 
previously deconvolved and beam-restored to the nominal resolution of our L'-band images). Color contours represent 
the bow-shock model projected onto the plane of the sky. PA = 55° (E of N) and 610=30° in all cases. We show the pure 
scattering case in orange (left), and equal contributions from scattering and thermal emission in purple (right). For each 
set of parameters, we show the results for two different values of i?o, in order to depict how change in the bow-shock size 
affects the fit. Contour levels are at (20, 30, 40, 50, 70, 90)% of the maximum. North is up, east is to the left. The axes 
are labeled in pixel offsets from the brightest point in the image; 1 pixel is equivalent to 0.027 mas. 



shaded area represents the disk of clockwise-rotating s t ars 



(CWS; IPaumard et al.||2006[ [Beloborodov et al. 



et al. 2009p , and the colored spheres are the ste 



2006 Lu 



ar mem- 

bers. The positions of the stars and the disk parameters 
{nx,n y, nz)=(-0.12, -0.79, 0.6) are from Paumard et al. 
(2006). The stars are represented by different colors ac- 
cording to their distance from the observer (green is closer 
and violet is further away from us). 



Eckart et al. (2002) show how a three-dimensional sep- 



aration r from the center can be estimated by comparing 
the proper motion (Vpm) of a star to the three-dimensional 
velocity dispersion a. The probability that a star at the 



position r has a proper motion larger than VpM can be 
calculated via 



PiV > VpM, 



= 1 - 



Jo 



V exp 



dv. 



(6) 



The velocity dispersion cr(r), and therefore also the proba- 
bility P(V>VpM,r), decrease with increasing radii r. For a 
fast-moving star it becomes increasingly unlikely that they 
belong to statistical samples at correspondingly larger radii. 
Therefore P(V>VpM,r) can be interpreted as a measure of 
how likely it is that the star belongs to a sample of stars at 
radius r or larger. Using the projected position R instead 
of the 3D value r, one can calculate P{V > VpM,R) as 
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an upper limit of the probability P{V > VpM,r). For a 
given radius r we can calculate the velocity dispersion as 
cr^ = GMbh/t, where Mbh is the mass of SgrA*. This 
gives a good estimate of a within the central parsec where 
the dynamics are dominated by the gravitational potential 
of SgrA*. Using this method, we estimate the line-of-sight 
position of X7 to be between z — —3.2" and z = 3.2", with 
P=67%. The elongation of the red object in Fig. [6] reflects 
this range of possible z-positions for X7. The 3D position of 
the bow-shock X3 cannot be constrained that well due to 
its lower proper motion and larger distance from the center. 
Therefore the orange object representing X3 is stretched 
across the entire z-axis. 

The mini-cavity is shown in pink: in projection we 
approximate it with a circle of ~1.2" radius. Paumard] 
let al.| (|2004) argue that the mini-cavity is a part of the 
NortEern Arm of the Mini-spiral. In their reconstruction 
of the streamer 3D morphology, the mini-cavity should be 
exactly in the plane parallel to the plane of the sky and 
containing SgrA*, or within ^2.5" from it. This is in ex- 



cellent agreement with orbital models of Zhao et al. (20091. 
The ±2.5" uncertainty of the line-of-sight position of the 
mini-cavity is depicted as an elongation of the pink object 
in Fig. |6] along the z-axis. 



7. Discussion 

All the interesting physics that can result from the model- 
ing is contained in the equation for the bow shock standoff 
distance Rq (eq.pl). The main drawback is a lack of knowl- 
edge about the stellar types (i.e. stellar wind parameters) of 
the two stars. The L'-band images are dominated by ther- 
mal emission of dust. Both X3 and X7 are bright at 3.8 fim, 
but extremely faint at shorter wavelengths. Both features 
are apparently not embedded in the mini-spiral material 
and therefore we expect them to be intrinsically dusty. 

If there was no influence of any external wind, the bow 
shock in our images would point in the direction of the 
proper motion. We understand that the Va in equations ^ 
and (|3| is the relative velocity between the stellar veloc- 
ity vector and the external wind velocity vector. Since we 
do not have full three-dimensional information about the 
stellar velocity, but only proper motions, the values of Va 
that we calculate in the following represent lower limits for 
the real external wind velocity. For the orientation of the 
projected vector Va we assume the position angle PA of the 
bow-shock symmetry axis (Figs. |4] and [5|. By estimating Va 
we can therefore determine the projected position of the ex- 
ternal wind source responsible for the observed bow shock 
morphology. 



7.1. Nature of the source at the position of X7 

In the following we discuss the possible nature of the 
X7 star, based on its apparent brightness in the K-band 
(16.9±0.1), combined with the distance modulus of 14.5 
and extinction of Ak «2.5. We discuss only stellar types 
that are in agreement with the K-band photometry. 

For the number density of the ambient medium we 



assume rt = 26cm~^ (Baganoff et al. 2003). This value 
is reasonable since X7 is not moving through the denser 
mini-spiral material. 



(1) Late B-type main sequence star (B7-8V). For O 
and B galactic stars the following relation holds: 

logim^v^Rl^) -1.37 + 2.07/og(L^/10^), (7) 



where rriu, is in M(;) yr 



in kms ^, i?* in Rq and 



in L0 (Lamers & Cassinelli 1999). For a typical 
B7V star with R* ~ 4Rq and LT^IO^ '^Lo we then 
have rn,„v,„ ~1O~®M0 yr""'^ kms~^. This kind of weak 

a negligible 



wind would require Va of only 15kms^^, 
velocity when compared to the proper motion of X7. Also, 
B-type main sequence stars are dust-free objects and thus 
probably not good candidates to produce features like X7. 



(2) Herbig Ae/Be (HAe/Be) star. This class of 
intermediate-mass pre-main sequence objects is charac- 
terized by strong wind activity and infrared excess. Line 
emission of HAe/Be stars often shows prominent P-Cygni 
profiles, indicating powerful winds with mass loss rates 
varying from 10~® to several times 1O~^M yr~^, and 
wind velocit i es of several hundred km s~^ (e.g. Benedettinil 

19951 



et al 
Nisini et 
between 



al 



1998 Bouret & Catala 1998 



|T995l 

loss 



r 



^ Nisini et al ^ 

show that there is a correlation 
rate and bolometric luminosity of 



mass 

HAe/Be stars. A star with L*=1O^'^L0 would have 
rriw ~ 1O~^M0 yr^^. Combined with ~500kms~^, 
this leads to Va ~3000kms~^. The external wind is 
then blowing with the velocity Vext ~2800kms^^ from 
the direction ^^60° (E of N). We have to note that the 
existence of the pre-main sequence stars at the GC is not 
established. It is not yet clear how the young (few Myr old) 
population in the central parsec has been formed, since 
the tidal field of Sgr A * prevents "norma l " star formation 



via cloud c ollapse (e.g. Nayakshin 2006b Portegies Zwart 
et al.||2006). However, the possibility of YSU presence in 



the central parsec was discussed by |Eckart et al.| (|2004|) 
andlMuzic et al.|(|2008|) 



(3) Central stars of planetary nebulae (CSPN). 

When all the hydrogen has been exhausted and the 
helium ignites in the core of a low/intermediate mass star 
(I-8M0), the star reaches the asymptotic giant branch 
(AGB). We note here that a typical AGB star at the 
GC is ~4 magnitudes brighter in the K-band than X7. 
AGB stars experience high mass-loss rates which remove 
most of the stellar envelope, leaving behind a stellar core. 
The star has now entered the post- AGB and subsequently 
planetary nebula (PN) phase. During this phase the central 
star remains at constant luminosity, while Teff rises, i.e. 
the star moves towards the left side of the HR diagram 
in a horizontal line. CSPN are characterized by small 
mass losses, but very fast stellar winds. A typical CSPN 



with m„ 



lO~'^M0yr-\ V . 



-3000kms-\ and R-R, 



■o 

(Lamers & Cassinelli 1999), will have L=3xlO'*L0 and 
therefore mif=17±l for logTe//~ (4.6-4.8). The evolution 
of a star from the AGB phase tow ards the final white dwarf 
stage is rapid (e.g. |Blocker| |1995 ), and the star will remain 
at the required brightness for less than 1000 yr before it 
becomes too faint. Adopting the typical CSPN parameters, 
we obtain Va ^2300 kms~^, and Vext ^2100 kms^^ from 
the direction -62° (E of N). 

(4) Low-luminosity Wolf-Rayet (WR) star; [WC]- 
type star. Population I WR stars are massive stars in a 
late stage of evolution. With typical luminosities between 



K. Muzic et. al.: Cometary shaped sources at the GC 



7 





4 2 -2 -4 
X (arcsec) 




^ 

a> 

en 
U 

^ 
>, 

s- 

a- 



2 0-2 

z (arcsec) 

Fig. 6. Three-dimensional view of some of the Galactic Center featm-es. The axes show offsets from Sgr A* (black sphere) 
in arcseconds. On the z-axis, positive means further away from the observer than Sgr A*. The shaded area represent the 
CWS disk and the colored spheres stars belonging to it; color scheme reflects the distance from the observer, with green 
being closest and violet farthest away from us. The bow shock sources are shown in red (XT) and orange (X3). Elongation 
along the z-axis reflects the uncertainty in the position of the two sources along the line of sight (see text). The pink 
spheroid represents the mini-cavity: in projection we plot it as a circle with radius r ^ 1.2" and the elong ation along the 
z-axis reflects the range of possible positions, as given by Paumard et al. (2004) and Zhao et al. (2009). We show the 
same setup from four different angles; the top right panel shows the projection onto the plane of the sky. 



3x10''' and 10^ Lq, they are several magnitudes brighter 
than X7, unless the local extinction at the position of X7 
is much higher than the estimated value of Ak=2.5. A 
class of CSPNe identified as WR stars share the origin with 
other PNe stars, but are spectroscopically classified as WR 
stars because emission-line spectra for both types of ob- 
jects come from expanding H-poor stellar winds. WR stars 
in P Ne almost exclusiv ely appear to be of the WC type 
(e.g. Gorny et al. 19951, and are labeled as [WC]. About 
15% of all the observed CSPNe appear to be [WC] ([Acker 



& Neiner 2003 ) . PNe stars with WR stars in their centers 
have, in gene ral, larger infrared excess th an normal CSPNe 



Crowther (2008) 



(!KwokIl2000!), and more powerful winds 
reviewed the properties of [WC] stars. Typically, winds are 
characterized by '-^ lO^^-lO^^M© yr^^ and ^-^200- 
2000 km s~^. Adopting the lower limit values for m^, and 
Vw, we obtain Va '^1900kms^^, Vext ^1700 kms~^ from 
the direction ^65° (E of N). For the average values m^, ~ 
lO"*"-^ and --1000kms"\ we have Va ~7700kms~\ 
Vext ^7400 km s^^ from the direction --55° (E of N). 
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7.2. Nature of the source at the position of X3 

In contrast to X7, there is no detectable point source at the 
position of X3 in our Kg-band images. Due to crowding, 
individual stars of Kg sal8 are difficult to detect within 
the central parsec and therefore we define this value as an 
upper limit on the brightness of X3 in the K-band. Here we 
discuss the possible nature of the star at the position of X3, 
and give possible solutions for the external wind direction 
and velocity. 

(1) Main sequence star. If on the main sequence, a star 
with K>18 would be of type AO or later. The mass-loss 
rate of low mass stars during the core H-burning phase is 



very small, of the order of 10~ - 10~ yr~ fLamers 
& Cassinelli 1999). In this case, any ISM streaming with 
the velocity of the order lOkms"^ could produce the bow 
shock with the standoff distance of X3. As in the case of 
X7, a main sequence star cannot be a source of a dust-rich 
envelope observed in the L'-band, and therefore we can rule 
out the main sequence nature of X3. 

(2) CSPN or [WC]-type star. As a CSPN evolves on 
the horizontal track in the HR diagram, it maintains the 
same luminosity while increasing the effective temperature. 
The peak emission of the stellar black body shifts blue- 
wards and the star becomes fainter in the infrared. A typ- 
ical CSPN or [WC]-star discussed above in case of X7, 
would have rriK >18 for logTej^>4.8. Adopting the typ- 
ical CSPN parameters, we obtain Va '-^ 1520 km s^"'^, and 
Vext ~1540knis^^ from the direction ^61° (E of N). In 
case of [WR]-star, using lower limit values as above, we 
obtain Va ^1240kms^^, and Vext ^1260 kms""'^ from the 
direction -62° (E of N). 

(3) Dust blob. Since with the current sensitivity of our 
observations we do not detect a point source in the K-band, 
we cannot rule out the possibility that X3 is a dust blob 
ablated by a wind from the direction of Sgr A*, rather than 
a stellar source. This scenario could explain the observed 
elongated tail. The proper motion suggests that the feature 
moves along with the rest of the material in the IRS13 
region. In this case modeling of the feature as bow shock is 
superfluous. 

7.3. Nature of the external wind 

The first question that we have to ask is about a wind 
source that can drive a shock of a certain velocity over a 
distance of few tenths of parsec. The argument of the ram 
pressure of such a wind at the given distance d leads to 



PaVs 



47rd2 



(8) 



Mw and are the mass loss-rate and velocity of the wind. 
The assumption behind Eq.[8] is that the wind emanates 
isotropically from a point source. Let us consider the case 
of a single star with the isotropic wind and typical wind 
parameters of the GC mass-losing stars: Miu=lO~^M0 yr~^ 
and Utu=750kms~^. This wind would be capable of driv- 
ing a shock of speed Vg into a medium of number den- 
sity uh with uhv'^ «1.7x IQ'^ x{d["])~^cm~^'km'^s~'^ , i.e. 
Ws=160 (800) kms~^ at a distance d=5" (1"), assuming 
n^f=26cm^'^. Note that this is an upper limit for Vg, since 
we have assumed that the ISM contains only atomic hy- 
drogen. These shock speeds are comparable to the proper 



motion values of the two sources and it is unlikely that 
the resulting bow shock orientations would be along the 
same line. Therefore it is clear that a single star cannot 
be responsible for X3 and X7. On the other hand, a com- 
bined isotropic wind of all the mass-losing stars at the 
GC (M^^lO^-'^Mgyr^i and i;^=750kms-i), would result 
in nnv'^g 1.7 x 10^ x{d["])-'^cin-^km'^s-'^ , i.e. 'i;s=1600 
(8000) kms^^ at a distance d=5" (1"). Interestingly, these 
results match the above estimate of X3 and X7 being 
CSPNe or [WC]-stars. Also, both types of stars are very 
short-lived, which nicely explains the lack of other objects 
of the same morphology in this region. 

The number of PNe expected to reside within the cen- 
tral parsec can be estimated as follows: 



PN 



E 

m<8M0 



iV*(m) 



TPN 

T{m) ' 



(9) 



where the summation goes over all stars with masses below 
8 M0 . We set the lower mass limit at about 1 M© , since 
the stars of lower masses have lifetimes comparable to, or 
longer than the Hubble time, and could not yet reach the 
PN phase. T(m) is a typical mass-dependent stellar life- 
time, and TpN is lifetime of the fast wind of the CSPNe. 
The latter is also mass-dependent, and ranges from <100 yr 
for a 8 MfT) star to ~ 3x 10^ yr for a 1 Mq (see e.g. Villaver 
et al.||2002 ). N^,{m) stands for the number of stars of the 
mass m currently present in the GC population. [Figer et al.| 
( 2004 ) present the K-band luminosity function (LF) result- 
ing from deep observations of the several HST/NICMOS 
fields in the GC region. To estimate the stellar counts at the 
faint end of the LF, we fit a power law to the number counts 
in the range Kwl6.5-19.5, where the data are reasonably 
complete. We also attempt to correct for the completeness 
using the data from Figer et al. ( 1999 1. By fitting the LF in 
this range we avoid the red clump population, but also the 
population of massive bright stars that dominate the LF in 
the central parsec. For scaling, we use the observed (com- 
pleteness c orrected) number co unts in the central parsec, 
at K = 17 (Schodel et al. 2007), while account i ng fo r the 
different magnitu de bin sizes in Schodel et al. (20071) and 
Figer et al] ( 2004[ ). This results in probability of ~40% to 
find one PNe within this region. We observe two sources in 
a much smaller volume, which makes it less likely that they 
are both of the same short-living type. This speaks in favor 
of X3 actually being a dust feature. 

The alignment of the two features is still not explained. 
In case of an isotropic wind arising from the mass-losing 
stars, one would expect a more random distribution of such 
sources around the center. Curiously, the two sources are 
arranged in the exact direction in which the mini-cavity is 
projected onto the sky. If we do not think of this arrange- 
ment as a chance configuration, this might indicate that (i) 
all three features (X3, X7 and the mini-cavity) are produced 
by the same event, and {ii) there is a preferential direction 
in which the mass is expelled at the GC. The possibility of 
a colli mated outflow was already discussed by |Muzic et al.| 
( 2007 ) . This outflow could also account for narrow dust fil- 



aments of the Northern Arm of the mini-spiral, as well as 
the H2-bright lobes of the circum- nuclear disc (CND). As 
the authors argue, the outflow could be linked to the plane 
of the mass-losing stars in the way that the matter provided 
by stars and not accreted onto Sgr A* is expelled perpen- 
dicular to the plane. Having an opening angle of about 30°, 
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this outflow could account for the mini-cavity, X3 and X7 
at the same time. In this case X3 and X7 should be lo- 
cated not too far away from the plane containing SgrA*, 
which is already suggested by the high inclination (iw90°) 
of the two bow shocks to the line of sight resulting from 
our modeling. 

8. Summary and conclusions 

We have presented L'-band observations of the two 
cometary-shaped sources in the vicinity of Sgr A*, named 
X3 and X7. The symmetry axes of the two sources are 
aligned within 5° in the plane of the sky and the tips of 
their bow-shocks point towards Sgr A*. Our measurements 
show that the proper motion vectors of both features are 
pointing in directions more than 45° away from the line that 
connects them with Sgr A*. Proper motion velocities are 
high, of the order of several 100 km s"""". This misalignment 
of the bow-shock symmetry axes and their proper motion 
vectors, together with high proper motions, suggest that 
the bow-shocks must be produced by an interaction with 
some external strong wind, possibly coming from Sgr A*, or 
stars in its vicinity. We have developed a bow-shock model 
in order to fit the observed morphology and constrain the 
source of the external wind. The stellar types of the two 
stars are not known. Moreover, one of the features is likely 
not a star, but just a dust structure. It might be located at 
the edge of the mini-cavity, and shaped by the same wind 
that produces the XT bow-shock. We discuss the nature of 
the external wind and show that neither one of the features 
can arise via interaction with an external wind originat- 
ing from a single, mass-losing star. Instead, the observed 
properties of the bow-shocks provide evidence for interac- 
tion with a fast and strong wind produced probably by an 
ensemble of mass-losing sources. Alternatively, a possible 
source of the wind could be Sgr A*. Shock velocities that 
can result from such a combined outflow over a distance 
assumed for the two features X3 and X7, match the ve- 
locities required to produce the bow shocks of stars in the 
late evolution stages of CSPNe or [WC]-stars. Short life- 
times of such stars can explain the lack of other similar 
cometary sources in the central parsec. We discuss our re- 
sults in the l ight of the par t ially-c ollimated outflow already 



proposed in Muzic et al. (2007) and argue that such an 
outflow, arising perpendicular to the CWS, can account for 
X3 and X7, as well as for the mini-cavity. The collective 
wind from the CWS has a scale of ~10 arcsec. On scales 
of about an arcsecond or less theoretical studies predict a 



radius-depen dent accretion flow (e.g. Narayan et al 
Yuan et al. [2003). Within this regiori the flo 



1998 



a major 



portion of tfie material originally bound for accretion onto 
Sgr A* is inverted and the material is expelled again to- 
wards larger radii. The presence of a strong outbound wind 
at projected distances from Sgr A* of only 0.8" (X7) with 
a mass load of IO^^Mq yr"-*- is in fact in agreement with 
models that predict a highly inefficient accretion onto the 
central BH due to a strongly radius dependent accretion 
flow. 

Knowledge of wind parameters for two bow-shock stars 
is crucial for drawing more quantitative conclusions on the 
nature of the external outflow. Spectroscopy should there- 
fore be the next step to confirm our hypothesis. 
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